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1. INTRODUCTION 

Wind energy (WE) is considered as the first sources exploited by humans after the energy of wood 
[1], [2]. It has been used to generate electrical energy for several decades [3], [4]. According to the renewable 
energies observatory, WE are at this time the most dynamic sector energy in the world [5]. Following the 
global oil crisis, the development and marketing of wind turbines was strongly encouraged [6]. Since then, 
the exploitation of wind resources has become increasingly efficient and the wind industry has experienced 
considerable growth over the last decade [7]. There exist two structures of wind turbines [3], [5]; the first 
structure operates at a fixed speed. It is directly connected to the electrical energy network. Whereas, the 
second structure operates at variable speed, this type of turbines can follow wind speed variations, such 
systems requiring power converters between generators and the grid. 

Several types of generators can be used such as asynchronous machine, dual feed induction machine 
and permanent magnet synchronous generator (PMSG). Newly, the exploitation of the PMSG becomes more 
and more common thanks to its advantages. Indeed, the direct connection of the PMSG to the turbine allows to 
obtain a significant at reduced speed [7]. For small power wind turbines, the use of a PMSG avoids the use of 
the gearbox, while for big power wind turbines, a gearbox is needed despite the use of the PMSG. The latter is 
connected to a rectifier, which supplies a DC bus connected on the other side to an inverter. This allows to 
separate the controls of the converters of stator and grid sides. For stator side converter (SSC), many control 
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methods are used in the literature, such as field-oriented control approach [8], feedback linearization control [9], 
these controllers are simple but they are classical and lack performance. Other works used sliding mode control 
[10] and backstepping control [11] however their main drawbacks are the complexity and difficulty of 
implementation. Thus, to overcome the drawbacks of the controllers mentioned above and thanks to the 
advantages given by Takagi Sugeno (TS) fuzzy controller applied to nonlinear systems and which has proven its 
performance and efficiency, we used in this work the TS controller and we applied the theory of Lyapunov 
stability. 

The aim of the control of grid side converter (GSC), makes it possible to control the production of 
active and reactive power as well as the maintenance of a constant voltage at the level of the DC bus [7]. In 
this part, we used the classical proportional integral (PI) controller thanks to its simplicity, so two control 
loops have been used in this section, an external control loop for the voltage of the DC bus, and an internal 
control loop for the powers. Thanks to the proportionality relationship between powers and currents, the 
power control is brought to a current control. So decoupling control of direct and quadrature currents is vital 
in order to control separetly the two types of powers [7]. 

The main contribution in this paper consists in using a TS fuzzy estimator and a TS fuzzy controller. 
The fuzzy estimator makes it possible, from any value of the wind speed, to predict the angular speed and the 
power converted at the level of the synchronous generator. While the TS fuzzy controller provides angular 
speed control, the electric power as well as the direct and quadrature currents of the synchronous generator. 
The paper is structured as shown in: The second section presents a general overview of the WECC and 
models each component. Section 3 describes the controller used for each part of the chain. The results of 
simulation are given in section 4 to confirm effectiveness of controllers. The paper is closed by a conclusion. 


2. WIND ENERGY CONVERSION CHAIN BASED ON PMSG 

The wind energy conversion chain (WECC) is made up of two parts, the first is defined by a turbine, 
a generator and a rectifier, while the second is made up of a DC bus, an inverter and a filter which is 
connected to the grid. The conversion principle described above is illustrated by the Figure 1. 
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Figure 1. General structure of wind energy conversion chain 


2.1. Wind turbine modeling 
The turbine is the first component of the WECC, the latter tranforms kinetic energy into mechanical 
energy. The converted aerodynamic power from wind is defined as [12]-[14]: 
1 
P= 5 ToMR*Cy(A)V™ (1) 
Where r,, C,, R, V and represent respectively the air density, the power coefficient, the length of the blade 
and the wind speed.The power coefficient is characterized by [15], [16]: 


—21 


Cy(A, B) = 0.5176 (116+ — 0.48 —5)e% + 0.00684 (2) 


Where 4; is described by as shown in [17] ,[18]: 
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A: 0.035 (3) 
Ai  A+0.08B B341 








A=— (4) 


Where, Q represents the mechanical rotation speed. The overall mechanical dynamics of WECC is given by 
the as shown in [15], [19]: 


J = Tn — Tem — fA (5) 


Where J represents the global inertia moment which depends on the turbine inertia, the generator inertia and 
the gearbox reduction. 


2.2. PMSG modeling 
The turbine is followed by a PMSG. From this machine, we obtain electrical energy. The 
electromagnetic torque of the machine in the d-q rotating frame is expressed as [20]: 


3 
Tem = 5 PLP els.q Piles Ls qs.als.q (6) 
For the PMSG we consider that [10]: 
Lsa = Ls.q = Ls (7) 


So, 
3 
Tem = 5 PEflsq (8) 


The dynamics behavior of the PMSG can be modeled by the as shown of state [21] : 
x(t) = Ax(t) + Bu(t) (9) 


Where A represents the state matrix, B the control vector, x describes the state vector and u is the control 
input. The state vector and the control input are respectively defined as: 


x(t) = [İsa isq Q] ,u= [Vsa Vogl" 


PQ 0 

ig 
_|_pg _® Poe 
A= PQ z rz 
9 er f 
2J J 


The state matrix A includes the angular velocity Q. However, the TS fuzzy model includes two local models 
each one is characterized by the state matrix A;. 








R R 
= Plan O — On © 
S S 
Rs POF Rs POF 
A, = =P nn ie P A A> = =P nax Ee A 
0 —3P QF _ f 0 —3P QF _ f 
2J J 2J J 


The TS fuzzy model of the system can be described by if-then rules. 


Rulei:If z(t)isF Then x(t)=A,x(t)+Bu(t) fori=1,2 (10) 
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The member ship functions are chosen as rectangular form, which nonlinearity is presented by two models: 


— N-Nmin 


U (11) 
F=1-K =p (12) 

The TS fuzzy model can be written as [22], [23]: 
X(t) = Dh mi(z@)) (Aix() + Bu(t)) (13) 


Where, r is the rule’s number (number of linear local model). The activation’s degree for rule is defined as: 
Fi 3 
ulz(t)) = = 0 <2) < 1,i=]1, r 
e Fj 


2.3. Grid side inverter modeling 
The two parts of WECS are linked through a continuous bus as shown by Figure 2. The modeling of 
the second part of the chain is given by these different equations [24]: 





_ i digd . 
Vga = Va — Rjiga — bp 24 + wgLpiga (14) 
_ digg 
V.a = Va — Retg.g — Lf Get 2gLfig.a (15) 
du 
C ma = lacs — lac.g (16) 


Where (Vsa V¢,) are grid voltages of d-q axes, (Vz, V,) designate the direct and quadrature components of 
the output inverter, (J, 4, Leq) represent the grid currents along d and q axes, Ry and L;represent respectively 
the resistance and the inductance of the filter, œ, represents the angular frequency, C is the capacitance of the 
capacitor, Uge represents the DC voltage of the capacitor, J,,.,1s the current of the SSC output to the DC link 
and lace is the current of the DC link output to the grid. The following two equations represent respectively 
the formulas relating to the active and reactive power [18]: 


P = Vo.alg.a + Vo.q lg.q (17) 


Q = Vo.qtg.a — Yg.atg.q (18) 


ides kiew 
Generator ee Grid 
E h ‘ 
Side ' | dc Side 
converter converter 


Figure 2. DC link capacitor 


3. CONTROL STRATEGIES OF THE WIND SYSTEM 


In this section, we are interested in controlling the WECS. It has two parts, the first relates to the 
control of the SSC and the second relates to the control of the GSC. 


3.1. Controller for generator side converter 

To maximize the conversion of wind energy, some equations of the system are used. The power 
coefficient 1s keeped at its maximum value Cpmax as well as the optimum value of the tip speed ratio. 
Likewise, mechanical rotational speed Q must be maintained at an optimum value and vary proportionally 
with the wind speed V. So, the mechanical reference speed is expressed as: 
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AoptV 
R 


(19) 





(re f= 
The power extracted by the wind turbine is also maintained at its maximum and it is expressed as shown in: 


t alre 
poo e (20) 
2o Àopt pmax 
max 
The direct reference current is equal to zero and the q-axis reference current is obtained based as shown in 


(10): 


_ 2Temref 
lsqref = 3 po, P (21) 


It is indispensable to control the mechanical speed of the turbine blades as well as the currents of the 
generator to increase the efficiency of the energy transformation. The diagram of the control strategy is 
presented in Figure 3, which includes two fuzzy blocks. The first block is a Takagi-Sugeno fuzzy estimator 
while the second is the Takagi Sugeno fuzzy controller. 





rectilpet 


: f> ~< [IS fuzzy T 
-i h - ELES] 
s SS csiimiter C} 
i l 4 sèf 
Figure 3. Stator side converter control 


3.1.1. T-S Fuzzy estimator 

Based on the minimum and maximum wind speeds (Vmin Vmax) allowed by the turbine as well as the 
corresponding optimum angular speed (Qop, Qop:2) to each of these wind speeds and the corresponding 
optimum power. It is then possible to predict, using a TS type fuzzy supervisor, the optimal angular speed 
Qop and the electrical power generated Pop: for any wind speed located in the (Vmin Vmax) range. Membership 
functions are shown in Figure 4. 


V V 


min max 


Figure 4. Membership function 


Nref E Qont = My Qopt1 + M2 Qont2 (22) 
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Pref z Popt = NıPoptı + N2Popt2 (23) 


Where, N} = 


; V 
y, min y, =—max 
P pt1andPopt2 denote the optimum power for the maximum and minimum wind speed respectively. 


3.1.2. TS type fuzzy controller 

To control speed and stator current, we propose parallel-distributed compensation (PDC) controller, 
which is based on TS fuzzy theory. The PDC controller has been proposed to stabilize the nonlinear system. 
The PDC controller rules can be written in the form: 
Rule: if zz) is Fz; and z,(t) is Fpi, So u(t)= - ki x(t) i=1,...,r 
Which is a state feedback controller to the consequence part. The TS fuzzy controller is expressed by: 


u(t) = — Vier H(z) Kix (t) (24) 

With k; represents the gain for local model for i=/,..,r . The following expression defined the tracking error: 
ey (t) = x(t) — Xref (t) (25) 

The derivate expression is expressed by: 
C= X= Meer (26) 


In a tracking problem, the expression of the control law (24) becomes as shown in [9]: 


u(t) = — Xi- Mi(z@)) Kie, (t) (27) 
So, 

ê, = Pf- mi(z(t)) (Aie, + Bu) (28) 

êr = Diy 4-1 Hi Mj (Aier — Bayer) = Lier Djn1 Hily (Ai — BK; ey (29) 
Theorem: 


The error of the overall closed loop system converge to zero, if there exists a matrix X definite 
positive which satisfies the following matrix inequalities: 


XA,’ + A,X — M;" BT — BM, < 0, i=1..r (30) 
XA; + A;X + XA," + A;X — M,"BT — M;' B" — BM, — BM, < 0, i<j=<r (31) 
Where M, = K,X and X = P“ 


Proof: 
The quadratic Lyapunov candidate function is chosen as follows: 


V(e,) = e1 Pe, (32) 
Where, e, is the error and P is the symmetric matrix definite positive. The derivate of the Lyapunov function is: 
V(e,) = é, Pe, +e, Pé, (33) 


Using as shown in (29), its follows that: 


V = (Sh Die Hi by (Ai — BK;)e,) Pe, + e," P(Yie1 Dias Hi uy (Ai — BK; )e,) (34) 
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V = Èi- Èi- Mi Hj (e, (A; E BK;) Pe, + e,"P(A;i a BK;)e,) (35) 
V = Ela Dye mim, (e,” ((4:" — K;7B")P + P(A; — BK;)) er) (36) 
V = Dies Xj=1 Mi Hy (e, (Ai P — Kj" BTP + PA; — PBK; )e,) (37) 


to prove that as shown in (37) is negative we transform the stability problem into LMI problem. 


V = Dim? (e," (AiP — Ki" BTP + PA; — PBK, ey) + Dina Dia H; U; (er? (Ar P+ 38) 
Aj" P + PA; + PA; — K;' BTP — K;' BTP — PBK; — PBK;)e,) 

It is very clear that as shown in (38) is a bilinear matrix inequality (BLMI). However, to transform 
as shown in (38) into an LMI, we must multiply the inequality in the left and the right by P”. Then we make 
the following change of variable X= P”. Then we obtain theas shown in: 


V = Vier mi? (e-7 (XA; — M;" BT + AiX — BM,)e,) + Dina Di H; 4, (e,7(XA;" + XA;" + (39) 
A,X + AjX — M;' BT — M;' B" — BM, — BM;)e,) 
with M; = K; X 

We remarque that as shown in (39) is the sum of two terms. To guarantee that this latter is negative. 
Each term must be negative. However, we get the following equations which verify the proposed theorem. 


XA,’ + A,X — Mi" BT — BM, < 0, i=1r (40) 

XA;' + A;X + XA;' + AjX — M;' B" — M; B" — BM, — BM; < 0, i<j<r (41) 

V(e,) = é,' Pe, + e, Pé, (42) 
SO, . 

V(e,) <0 (43) 


3.2. Controller of grid side converter 

In this part, we are interested to control the DC voltage across DC bus and powers, which will be 
transmetted to power grid. According to as shown (17) and (18), the two powers previously defined are 
strongly coupled, which makes it difficult to control them independently. However, it is necessary to 
decouple them. For this, we choose to align the direct axis component of the grid voltage with the direct axis 
of the field-oriented control frame and we keep zeroing the componen of grid voltage along the quadrature 
axis [25]. 


Voq = 0 (44) 
Via = V, (45) 


Therefore, the powers equations become as shown in [25]: 

P = Viga (46) 

Q = Volga (47) 
Therefore, the expressions of these two powers are decoupled as shown (46) and (47). The currents along the 


direct and quadrature axes allow to controlrespectively the active and reactive powers. Figure 5 describes the 
diagram structure of the current controllers along the d and q axes. 
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Figure 5. (a) currents control along d-axis, and (b) currents control along q-axis 


The transfert function of the open loop is given by : 


__ KitKpS 1 
Hor = S (LfS+Rf) (48) 








To determine the parameters of the PI controller, we use the pole compensation method. However, we choose 
So, the transfer function in closed loop becomes as: 


1 


Ac, = ey (49) 
where, T = “L = 0.05 


We propose a PI type controller to keep a constant DC bus. In this case, two control loops are 
obtained along the direct axis. A first internal loop for controlling the current and a second external loop for 
the voltage of DC bus. In this case, the reference current along the direct axis is provided by the output of the 
controller. Figure 6 describes the general structure of the DC bus voltage controller. 





ee _——. 


Figure 6. General structure for DC bus control 


The current transfer function represents a first order function, where T << L. However, we 
approximate the function to a unit gain, so the transfer function of open loop of the DC bus voltage can be 
written as: 


Kj, 1 


The transfer function of closed loop is presented by: 


Nee fea. (51) 


Fi E 


This function is of the form: 


1 
H = —— 32 

—s2 428 cay ( ) 

WO WO 
With ¢ describes the damping factor and wo represents the own pulsation of the system. In this work ¢=/ and 
@o=40 rad/s. Finally, we obtain the parameters of the DC bus controller: 
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K; = Ca? (53) 
K, =2¢Ca@ (54) 


Figure 7 summarizes the GSC control part, indeed, the measured currents are compared to their 
references. The direct reference current is obtained from the DC bus controller, while the quadrature 
reference current is zero to obtain zero reactive power. Two PI controllers are proposed, one for direct 
current and the other for quadrature current. The compensation terms are added to the voltages obtained at 
the output of these two controllers. The two finally obtained voltages are transformed by the PARK 
transformation to obtain the modulating voltages which will be compared with a carrier: this is the principle 
of pulse width modulation (PWM). 


DC bus Inverter 


j Filter 
U, t= = RI 






[Park |B | Park | 
rT iy 


Pi 


Figure 7. Grid side inverter control 


4. RESULTS OF SIMULATION 
In order to chek the performance of the controllers, some results obtained are presented in this 


section. For this, the system parameters used are given in Table 1, Table 2 and Table 3. 


Table 1. Turbine parameters 
Radius of blades Tip speed ratio Pitch Angle Rated power 
R=3m Aop=8.1 B=0 P=10 Kw 


Table 2. PMSG parameters 


Inertia Friction coefficient Stator resistance d-q axis inductance Magnetic flux Number of poles pair Gearbox 
J=0.05N.m F=0.006 Nm sec/rad R, = 0.5 Ohm L,= 1 Mh D, = 0.16 wb P=8 G=8 


Table 3. DC bus parameters 
Capacitor Udc_ref 
10e-4 F 400 V 


The resolution of the LMIs (30) and (31) gave the following matrix X and gains K;;: 
X= [0.3230 -0.0000 0.0000; -0.0000 0.3284 -0.6670; 0.0000 -0.6670 58.4888] 
K,= [1.1295 -67.3718 136.8268; 66.2580 1.0497 5.7163] 
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K= [1.1295 -111.1048 225.6451; 109.2680 1.0497 5.7163] 
K; = [1.1295 -133.3102 270.7425; 131.1063 1.0497 5.7163] 
K4 = [1.1295 0.0089 -0.0180; -0.0087 1.0497 5.7163] 

The parameters of the PI controllers have been chosen as: 
For the currents: K,=0.05 and K;=10 

For the DC bus voltage: K,=0.8 and K;=16 

The real wind speed can take any value, for this work we chose as shown in curve Figures 8 tol0. 
The power coefficient reaches its maximum and it is maintained at this value despite the change of wind 
speed as given in the Figure 9. This coefficient guarantees the extraction of the optumum power as described 
in Figure 10, in fact the power is maintained at the maximum value for the two values of the wind thanks to 
the best value of power coefficient. 

The angular speed is given by Figure 11. It is clear that the measured angular speed follows the 
reference angular speed even when changing the wind speed. The evolution in the behavior of the currents 
along the axes q and d is presented in Figures 12 and 13. It is clear that the curves obtained by the controller 
follow the reference curves. In fact, the current obtained along the q-axis is near to its reference in the two 
wind intervals. Also, the current along the d-axis is near to zero in both intervals. 

Figure 14 demonstrates that the voltage across the DC bus is continuous and practically constant 
despite the change of the wind speed. The voltage value follows its reference. Figures 15 and 16 describe the 
obtained currents along axes d-q of the power grid side, these currents are near to the reference ones. The d- 
axis current follows its reference with a small overshoot when changing the wind speed. Also, the current 
along q-axis converge to zero. 

Figure 17 represents the currents obtained at the output of the filter. These currents are three-phase, 
balanced and offset by 2*z/3. Their amplitude is influenced by the change of the wind speed, in fact they 
have the same appearance, and this is clarified in Figure 18 which we zoomed in for the two wind intervals to 
show the amplitude of the current in these two intervals. 

Figure 19 represents the voltages obtained at the output of the filter. These voltages are three phases, 
balanced and offset by 2*z/3 and their amplitude is constant despite the change of the wind speed, this is 
clarified in Figure 20. Figure 19 represents the voltages obtained at the output of the filter. These voltages are 
three phases, balanced and offset by 2*z/3 and their amplitude is constant despite the change of the wind 
speed, this is clarified in Figure 20. 

Figures 21 and 22 demonstrate that the two powers previously defined follow their references. In 
fact, the active power is near to its reference with a fairly small delay while the reactive power is practically 
zero. The tracking of powers to their references is due to the proportionality relationship between the currents 
and the powers and confirms the performance of the controller applied for control of the currents. 
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Figure 18. (a) zoom on currents for the first interval; (b) zoom on currents for the second interval 
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Figure 19. 3-phase voltage at the filter output 
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Figure 20. (a) zoom on the voltages for the first interval; (b) zoom on the voltages for the second interval 
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Figure 21. Active power Figure 22. Reactive power 


5. CONCLUSION 

In this paper, we are concerned into control the WECC. However, we started with the modeling of 
the system. The latter has two parts, a first is called SSC and a second is called GSC. In the first part we have 
used a fuzzy estimator and a Takagi Sygeno type controller to maximize the conversion of wind energy. The 
fuzzy estimator makes it possible to predict the optimum mechanical rotation speed of the blades as well as 
the optimum power that the turbine can transform from the wind. In the second part, we have used a PI type 
controller. The goal is to retain the DC bus constant as well as the control powers. The results of simulation 
are given to demonstrate the performance of the proposed controllers in the two parts. It can be concluded 
that the two controllers have succeeded in achieving the desired objectives. 
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